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Abstract: We have performed molecular dynamics simulations of Fe§oNi2o alloys using 
an inter-atomic potential of the EAM-type which allows the simulation of the martensite- 
austenite transition. We present results, showing the development of an inhomogeneous 
shear system on a nanoscale during the thermally induced austenitic transformation. In 
addition to this we obtained the phonon dispersion relations of the martensite phase by 
calculating the dynamical structure factor from our simulation results. On approaching the 
transition temperature the phonon dispersion shows anomalies which might be connected 
with the formation of the microstructure during the austenitic transition. 



1. INTRODUCTION 

Relationships between atomistic processes and the formation of the microstructure during martensitic 
transformations are not well understood today. In order to get more insight into this we have 
performed molecular dynamics simulations of Fe x Nii_ x alloys. Within the range 0.66 < x < 1.0, 
these alloys show experimentally a martensitic transformation from a fee high-temperature phase to 
a low-temperature bec phase (see for example p| and references therein). 

Most investigations in the context of martensitic transformations concentrate on the austenite to 
martensite transition. This is reasonable since the formation of microstructure during the martensitic 
transformation influences the properties of the low-temperature phase to a high degree and makes 
it rather difficult to investigate the features of the homogeneous martensitic phase. Nevertheless 
it is desirable to study the homogeneous phase, in order to be able to distinguish between genuine 
features of the martensitic phase and those caused by the microstructure. Therefore, we have done 
simulations in the low-temperature bec phase, looking for the processes leading to the austenitic 
transition which also lead to the formation of a microstructure in the high-temperature phase. 

We present results of two distinct simulation sequences. We begin with a study of the formation 
of a microstructure on a nanometer lengthscale. These simulations required the consideration of a 
large number of atoms. In order to do this we employed a semi-empirical potential which is based 
on the embedded-atom method (EAM) introduced by Daw and Baskes ]2|, [|. This enabled us to 
simulate systems with characteristic linear dimension larger than 5 nm. Despite of the fact that 
this is, physically spoken, still a small system, we find the formation of a microstructure with an 
inhomogeneous shear system. 

In a second investigation we calculated the phonon dispersion curves of FesoNi 2 o in the bec 
phase. Our results show that the austenitic transition of FesoNi 20 is accompanied by similar phonon 



2. COMPUTATIONAL METHODS 



Classical molecular dynamics simulations have been performed, using a semi-empirical potential 
based on the EAM which has been constructed recently for the study of martensitic transformations 
in Fe x Ni!_ x alloys 0, |6|. In contrast to simple pair-potentials, potentials based on the EAM are able 
to describe the elastic behaviour of a metal correctly, while remaining computationally efficient. 

Standard techniques of molecular dynamics simulations H] like periodic boundary conditions and 
the Verlet algorithm for the integrations of the equations of motion with a time step 5t = 1.5 x 10~ 15 s 
have been employed. Both simulation sequences started with configurations of atoms on an ideal bcc 
lattice with a random distribution of Fe and Ni atoms. 

The first simulation sequence was done with a configuration of 16000 atoms (20 x 20 x 20 cubic 
bcc elementary cells) having an iron concentration x = 80.01 %. This system was heated from a 
temperature of T = 600 K to 700 K in steps of 50 K. Afterwards the temperature was increased in 
steps of 20 K until the austenitic transformation occurred at T = 860 K. At each temperature 1000 
simulation steps were done in order to reach thermal equilibrium and another 10000 steps (15 ps) for 
measurement purposes. These simulations were carried out within the isothermal-isobaric ensemble 



generated by the Nose-Hoover thermostat [|| § and the Parrinello-Rahman scheme |10| , |TT[ with a 
fluctuating simulation box. 

The second set of simulations used a smaller configuration of 12 x 12 x 12 cubic elementary 
cells (3456 atoms) with the same iron concentration x = 80.01%. This system was simulated at 
temperatures T = 300, 500 and 700 K with a fixed simulation box and within the microcanonical 
ensemble. From a previous investigation the equilibrium lattice parameters were known. After the 
equilibration phase 40000 simulation steps (60 ps) were done, writing the atomic positions to a file 



after every 10th simulation step. From this file the dynamic structure factor [12 



1 r dt 

S(q^) = jj E J -e^(exp[zq. U (R',0)]ex P Hq-u(R,t)]) (1) 

has been calculated (N is the number of atoms and u(R,£) represents the displacement of the atom 
from the ideal lattice position R at time t) . The phonon dispersion curves were obtained by evaluating 
the positions of the peaks of S(q,u) at different wave- vectors q. 



3. RESULTS 



3.1 Microstructure formation 

The simulations done with 16000 atoms revealed the same general characteristics of the martensitic 
phase and the austenitic transition as similar calculations of smaller systems did0, |(J. In particular 
identical structural and orientational relationships were observed. One of the (110)bcc plane sets 
changed to a set of close packed planes, remaining the [001] direction parallel to this planes unrotated. 
However, more interesting is the structure of the resulting austenite phase. Figure 1 shows two atomic 
layers of the resulting structure viewed along the unrotated [001] direction. The stacking sequence of 
the close packed planes turns out to be ABCABCA CABABCABCB CAC (from left to right). Keeping 
in mind the periodic boundary conditions, it can be seen from the stacking sequence and Fig. 1 that 
the whole system consists of two big fee plates (6 and 8 layers thick) which are separated by thin 3 
layer plates. The slip faults which separate these basic blocks lead to an inhomogeneous shear angle 
of 2.5°. The atomic movements that generated the structure of the system in Fig. 1 are elucidated by 
Fig. 2, which displays the deviations of the atoms drawn in Fig. 1 from ideal body-centered lattice 



Figure 1: Two atomic [001] layers of the simulated FesoNi 2 o crystal after the austenitic transformation. The diagonal 
lines indicate stacking faults of the (111) planes. 
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3.2 Phonon dispersion relations 

In order to see the effects of alloying on the phonon dispersion curves, first the phonon spectra of pure 
Fe and Ni were calculated analytically from the EAM-potentials by diagonalization of the dynamical 
matrix. Figure 3 shows the result, obtained for a bcc structure with lattice constant = 5.4369 Ob 
determined for FegoNi 20 at T = 300 K from the simulations. Along the q = [111] and [001] directions 
both elements show a normal behaviour. Along the other directions considered in Fig. 3 this is only 
true for Fe. Around q = [^ |l] and along the q = [110] direction the spectrum of Ni exhibits a strange 
behaviour and one of the modes has negative squares of frequencies here. This result itself is not 
very exceptional since experimentally Ni shows no stable bcc phase. But the unstable mode along 
the [110] direction corresponds in the limit q — > to a negative elastic constant C = |(Cn — Cu). 
This is consistent with ab initio electronic structure calculations which also reveal a negative value 
ofC'p. 

With respect to this instability of Ni and the fact that the bcc - fee transition requires a shear 
along [110]bccj this direction appears to be a good candidate for phonon anomalies in those Fe x Ni]_ x 
alloys, which exhibit a martensitic transformation. Figure 4 shows the phonon dispersion curves of 
Fe8oNi2o along this direction determined from molecular dynamics simulations by calculation of the 
dynamic structure factor S(q,u) at T = 300 K (a), and the temperature dependence of the low 
lying TA2 mode (b). The data points in these figures represent the positions of the peaks in S(q,u) 
averaged over all crystallographically equivalent [110] directions. The error bars are derived from the 
corresponding standard deviations. If no error bars are drawn, the standard deviations are smaller 
than the symbol sizes. 

Figure 4(a) demonstrates the dramatic effect of Ni on the phonon dispersion curves. The Ion- 




to this in Fig. 4(b) it can be seen without doubts that the low lying transversal mode develops two 
distinct anomalies on approaching the transition temperature. At q — > the slope of the branch 
reduces with increasing temperature, leading to a positive curvature of the corresponding dispersion 
relation. The second anomaly is located between q = 5 [110] and the Brillouin-zone boundary. The 
dip occurring here in the transversal mode is also accompanied by large errors of the data points. 



4. DISCUSSION AND CONCLUSIONS 



Results of molecular dynamics simulations shown in Fig. 1 and Fig. 2 demonstrate the formation of 
an inhomogeneous shear system during the austenitic bcc - fee transition in Fe x Nii_ x alloys. This 
shear system is quite similar to that observed in the austenite martensite - transitions of these 
systems. The microstructure consists of homogeneous fee plates separated by slip faults. Though the 
length scale of these simulations is still far too small to give definite answers, it is interesting to see 
that the thicker fee plates are separated by rather thin plates consisting of only three atomic layers. 
This could be just a random result of this particular simulation. But this observation is confirmed 
by the anomalies of the phonon dispersion curves resulting from the second simulation sequence, 
since wavelength and direction of the anomaly between q = |[110] and the zone boundary are in 
accordance with the microstructure observed in Fig. 1. Therefore we think that the formation of 
the microstructure during the transformation can be attributed to this anomaly which is visible at 
temperatures far below the actual transformation. Nevertheless, simulations of even larger systems 



From Fig. 4(a) the effects of the addition of Nickel on the phonon dispersion curves can be seen. The 
most striking effect is the reduction of the low-temperature value of C . This explains the decrease 
of the austenitic transition temperatures with increasing Ni contents [[!]. 

It is clear that under normal experimental conditions, the properties of the martensitic phase 
are dominated by the microstructure which develops during the martensitic transformation. But 
we think that the results of our present molecular dynamics simulations show that the austenitic 
transition also has its own genuine features which are worth to be studied. To our knowledge there 
are no experimental data of the phonon dispersion relations of the bcc phase of Fe x Ni!_ x available 
at the time being. It will be very interesting to see if our theoretical prediction can be confirmed 
experimentally. 
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